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COLLISIONLESS  EFFECTTS  ON  THE  SPECTRUM  OF  SECONDARY 
AURORAL  ELECTRONS  AT  LOW  ALTITUDES 


fi 


Introduction 

The  differential  flux  of  primary  and  secondary  electrons  under 
auroral  activity  has  been  the  subject  of  a series  of  experimental 
measurements  [Feldman  et  al.  1971:  Amoldy  and  Choy  1973:  Reasoner 
and  Chappel  1973;  Feldman  and  Doering  197^:  Matthews  et  al  1976].  A 
conmon  feature  of  all  these  rocket  experiments  is  that  between 
30-33  eV,  where  the  measurements  overlap  in  energy,  the  electron  flux 
data  can  be  fit  by  a power  law  E with  .3  ^ or  ^ 1.0.  Theoretical 
models  based  on  ionizing  collisions  of  the  primary  (2-30  kev)  electrons 
producing  secondary  and  backscattered  electrons,  predict  much  steeper 
( ~ E~^)  and  strongly  altitude  dependent  spectra  in  sharp  contrast 
with  the  observations  [Banks  et  al  197^:  Berger  et  al  1974;  Rees 
and  Maeda  1973] * It  should  be  noted  that  the  above  models  are  in 
reasonable  agreement  with  observed  fluxes  for  energies  below  20  eV  (Flg.l). 

In  an  earlier  paper  Papadopoulos  and  Coffey  (1974a)  demonstrated 

that  collisionless  processes  can  be  very  Important  in  the  auroral 

arcs  and  could  explain  the  characteristic  features  of  the  secondaries 

observed  by  Reasoner  and  Chappel  [l973]  At  altitudes  6OO-7OO  km.  More 

recently  Matthews  et  al  [ 1976]  demonstrated  that  similar  processes  can 

operate  under  Intense  auroral  conditions  down  to  altitudes  of  I30  km, 

and  could  account  for  several  of  their  experimental  observations. 

The  fundamental  physical  idea  of  the  Papadopoulos-Coffey  theory  was 

that  the  primary  electron  spectrum  has  the  form  of  a large  flux  spike. 

Note:  Manuscript  submitted  October  6,  1977. 
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which  can  thus  be  unstable  and  produce  plastoa  waves  with  phase  velocities 
near  the  velocity  o£  the  precipitating  particles.  When  these  waves 

reach  a certain  amplitude  they  become  themselves  unstable  to  parametric 
instabilities  and  decay  to  plasma  waves  of  low  phase  velocity  which 

can  interact  with  the  ambient  plasma  and  create  fast  ( > 20  ev) 
upstreaming  and  downstreaming  secondary  electron  fluxes.  An  additional 
effect  of  the  above  process  Is  the  excitation  of  low  frequency  Ion 
density  fluctuations  which  can  produce  anomalous  resistivity  in  the 
auroras  [Papadopoulos  and  Coffey  197^1>]  • Detailed  accounts  of  the 
theory  can  be  found  In  the  forementloned  papers  as  well  as  in 
Papadopoulos  C 19751  and  Papadopoulos  et  al  [ 197^1 • Th®  derivations 
presented  In  the  above  papers  are  totally  analytic  and  therefore 
several  approximations  were  employed.  In  view  of  the  Importance  of 
the  problem  In  the  understanding  of  the  auroral  processes,  we  decided 
to  test  some  aspects  of  the  formation  of  secondaries  by  using  particle 
computer  simulation  techniques,  which  allow  us  to  solve  the  complete 
non-linear  problem,  and  compare  the  results  with  the  experimental 
observations  of  Feldman  and  Ooerlng  [ 19751 • ^®  complete  simulation 

of  the  problem  should  Include  the  colllslonal  processes  via  a Fokker- 
Planck  collision  Integral  that  computes  energy  deposition  rates  In  a 
similar  fashion  as  In  Banks  et  al  [197^1  or  Walt  et  al  C1969I  and  the 
wave-particle  Interactions.  In  view  of  the  large  disparity  In  the 
time  scales  of  these  processes,  the  only  reasonable  way  of  producing 
such  a program  Is  by  describing  the  wave-particle  Interactions  by  a 
diffusion  tensor  in  velocity  space  whose  properties  have  been  determined 
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by  simulation  studies  of  the  collisionless  microphysics.  We  are  pre- 
sently developing  such  a capability  and  the  present  study  represents 
a step  in  this  direction. 
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In  this  paper  we  assume  that  the  spectrum  of  the  secondary 
electrons  is  of  the  form  E”®  (actually  e”®*^  to  simulate  the  Feldman- 
Doering  results)  as  expected  from  the  collisional  theory  and  examine 
the  modifications  that  will  be  introduced  due  to  the  wave-particle 
Interactions . 

The  present  study  starts  with  two  basic  assumptions: 

(a)  A secondairy  flux  of  the  form  e"®*^  has  been  established  by 
collisional  processes. 

(b)  Electron  plasma  waves  with  phase  velocities  near  the  velocity 
of  the  primary  spike  like  electron  spectrum  can  be  generated. 
Generation  of  such  waves  is  a basic  concept  in  plasma  f>hysics  due 
to  the  Interaction  of  precipitating  beams  with  the  plasma.  Such 
waves  have  been  recently  observed  in  the  aurora  at  altitudes  of 
180  km  and  with  a modulation  amplitude  5^  [Spiger  et  al  19751 • 

Using  particle  simulations,  where  10^  particles  are  followed  in  their 
self  consistent  orbits  and  where  the  above  assumptions  are  used  as 
inputs  at  the  initial  time,  we  try  to  answer  the  following  questions: 

(1)  Does  the  presence  of  the  plasma  waves,  which  are  non-resonant 
with  ambient  particles,  change  the  observed  flux  power  law? 

(ii)  At  what  energies  the  break  occurs? 

(iii) ls  this  energy  a function  of  the  wave  energy  level? 
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[iv)  What  is  the  final  high  frequency  and  low  frequency  wave 
energy  spectrum? 


r 


In  order  to  keep  as  close  as  possible  with  the  experimental  data  we 
carried  our  simulation  for  the  parameters  of  the  experiment  of  Feldman 
and  Doering  [197^1,  however  the  results  have  a much  broader  applicability. 

In  the  next  section  we  describe  briefly  the  numerical  code  used 
for  the  studies.  In  Section  III  we  present  the  results.  In  the  final 
section  we  discuss  the  constraints,  conclusions  and  the  Implications 
of  the  study. 

II.  Description  of  the  Code 

The  simulations  were  done  with  a one-dlmenslonal  pure  electro- 
static hybrid  code.  The  electrons  are  treated  as  discrete  particles, 
the  Ions  as  a fluid.  For  this  problem,  the  Ion  dynamics  can  be 
accurately  described  by  the  fluid  equations  [DuBoIs  and  Goldman  I967] « 

By  using  such  a representation  for  the  Ions,  we  reduced  the 
statistical  noise  in  the  Ion  background  due  to  the  finite  number  of 
simulation  particles.  By  reducing  the  statistical  fluctuation  level, 
we  were  able  to  simulate  much  weaker  Instabilities  then  would  have 
been  possible  with  a pure  particle  simulation.  Since  we  are  Inter- 
ested In  wave-particle  Interaction  Involving  the  electrons,  we  kept 
a particle  representation  of  the  electrons.  The  particle  side  of 
the  code  Is  based  upon  PPOWER  [Boris,  et  al.  1973]  ®nd  is  typical 
of  particle  simulation  codes  presently  In  use.  The  particles  are 
advanced  using  the  Internally  (via  Poisson's  eq.)  calculated  electric 
fields.  The  fluid  quanCltles-denslty,  momentum  and  energy  are 
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advanced  In  time.  A new  electric  field  is  then  calculated  from  the 
sum  of  the  fluid  and  particle  charge  density  and  the  process  is 
repeated.  The  fluid  algorithm  used  in  the  minimum  diffusion  SHASTAX 
[Book  et  al.  1975] • 

III.  Simulation  Results 

The  initial  configuration  consists  of  a Maxwellian  electron 
plasma  with  thermal  velocity  an  initial  symmetric  secondary  electron 
flux  which  starts  at  ~ 5^^  and  decreases  as  till  | v^^^j “ 6v^ 

with  a total  of  2.256  of  the  total  number  of  electrons.  The  ions  are 
treated  as  a fluid  with  a temperature  equal  to  the  electrons.  The 
plasma  waves  were  excited  with  phase  velocity  v^^  = 8l  V^.  In  actual 
numbers  assuming  ambient  electron  temperatures  of  the  order  of 
•5  eV,  the  initial  secondary  flux  extends  from  4.5  - 18  eV  and  the 
waves  are  created  by  a flux  with  a spike  at  energies  larger  than 
3.2  keV.  The  waves  were  allowed  to  grow  up  to  a certain  wave  energy 
level  Wi.  The  ratio  of  wave  energy  level  over  the  electron  thermal 
energy  Wi/nT^  at  which  simulations  were  performed  varied  between 
.1  - .5.  We  will  comment  on  the  significance  of  such  values  in  the 
next  section.  The  physical  basis  for  a model  where  the  beam  waves 
grow  up  to  a certain  level  and  are  subsequently  depleted  has  been 
described  in  Papadopoulos  [1975]  subsequently  confirmed  numerically 
in  Falmadesso  et  al.  [ 19751  Rowland  and  Papadopoulos  [ 19771  • 

The  results  are  shown  in  Figs.  1-4.  In  Fig.  1 we  show  the 
initial  and  final  electron  flux  for  the  three  cases  under  study.  As 
seen  from  Fig. 2 and  from  Table  I the  effect  of  the  collisionless 
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interactions  is  to  change  the  power  law  of  the  secondaries  from  E 


-3.2 


to  E - . The  value  at  which  the  break  occurs  depends  on  the  value  of 
U/nTg.  The  dependence  seems  to  be  in  accordance  with  the  theory 

[Papadopoulos  1975^  Matthews  et  al.  197^  which  predicts  that  the 
velocity  v^  at  which  the  break  occurs  scales  with  \ 


(Fig.  3)-  Unforttmately  due  to  computer  limitations  we  could  not 
perform  simulation  at  much  lower  values  of  W^nT^.  As  will  be  discussed 


in  the  next  section, the  values  of  used  correspond  to  strong 

precipitation  with  number  densities  of  the  order  of  1 cm~^.  Figure  k 
shows  the  energy  spectrum  of  the  electron  plasma  waves  as  function 
of  the  phase  velocity  at  the  initial  and  the  final  time.  It  is  seen 
that  the  waves  are  parametrically  coupled  Co  lower  phase  velocities 
where  they  interact  with  Che  ambient  particles  and  modify  their  dis- 
tribution function.  On  Che  basis  of  Che  observed  spectrum  one  can 
determine  the  diffusion  coefficient  in  velocity  space  which  will  be 
used  in  the  study  which  includes  both  collisional  and  collisionless 
effects.  Figure  5 shows  the  spectrum  of  the  low  frequency  density 
fluctuations.  From  this  spectrum  one  can  compute  the 
anomalous  resistivity  as  discussed  in  Papadopoulos  and  Coffey  [ 197^h] . 
IV.  Conclusions  and  Discussion 


On  Che  basis  of  Che  above  results  we  can  answer  the  questions 
set  forth  in  the  Introduction.  It  has  been  shown  that  Che  existence 
of  plasma  waves  even  in  a region  where  they  are  non-resonant  with  the 
ambient  particles  can  significantly  modify  the  observed  flux  power 
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law.  It  was  also  shown  that  while  the  index  of  the  power  law  is 


independent  of  the  initial  wave  energy  density,  the  velocity  at  which 

Vfa  j nf 

the  break  occurs  is  consistent  with  the  relationship  — — . 

V 1 Wi 
e 

Finally  large  density  fluctuations  were  created  by  the  ponderomotive 
force  of  the  high  frequency  waves.  Such  waves  can  produce  anomalous 
resistivity  as  discussed  in  Papadopoulos  and  Coffey  [197^]  and  Rowland 
et  al.  [1977]. 

Applicability  of  the  Model  to  the  Aurora  Zone 

While  the  model  as  discussed  above  and  therefore  its  conclusions, 

is  completely  self  consistent,  its  applicability  to  the  auroral  zones 

requires  some  further  discussion.  What  we  have  demonstrated  here  is 

that  given  an  initial  secondary  spectrum  such  as  expected  from 

Wi 

collisional  interactions  and  a wave  energy  level  — =r-  .1  - .5/  one 

^^e 

can  account  for  the  observations.  Electron  plasma  waves  have  been 
detected  in  the  auroral  zones  with  energy  densities  comparable  to  the 
above  [Spiger  et  al  19751  • addition  observation  of  the  primary 
electron  spectrum  indicates  that  such  waves  must  be  generated  in 
accordance  with  plasma  theory  [ Papadopoulos  and  Coffey  197^a»  Matthews 
et  al  197^ • As  discussed  in  the  previous  papers  if 

where  e the  energy  of  the  precipitating  particles  and  V , AV  their 

velocity  and  velocity  spread  the  value  of  — will  be  given  by 
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Wi 


nT 


» 10^  A 


(2) 


where  A is  the  atomic  number  of  the  dominant  ionic  species.  Inequality 
(1)  is  easily  satisfied  in  the  auroral  zones.  It  is  important  to 

Wi 


notice  that  in  £q.  (2) 


nT 


does  not  depend  on  the  energy  of  the 


precipitating  particles  but  only  on  the  relative  energy  spread  which 

observations  indicate  is  relatively  constant.  Taking  A **  l6, 
v_ 

n *=«  10®cm  ® and 
Wi 

' 10‘^"  ^ 


AV 


5 as  rather  constant  conditions,  it  turns  out 


that 


nT 


n^^^.  The  value  of  n^  is  the  only  one  that  can  vary 


significantly.  For  strong  auroras  n^  1 and  the  values  used  in  our 
simulations  are  appropriate.  For  weaker  auroras  n^  ~ 10“^  and  much 
smaller  values  of  should  be  used.  Unfortunately  the  simulations  for 
such  values  become  increasingly  lengthy  and  we  could  not  carry  them 
out.  However  we  expect  the  basic  results  to  be  valid  since  even  for 
such  weak  fields  the  conditions  for  parametric  coupling  to  secondary 
waves  are  easily  satisfied. 

In  conclusion  we  can  state  that  the  particle  simulations  confirm 

the  fact  that  collisionless  interactions  can  be  responsible  for  the 

break  from  e”®  dependence  of  the  secondary  electrons.  The  velocity 

Vb  I nT" 

at  which  the  break  occurs  is  given  by  — . For  the 

above  parameters  we  find  the  scaling  law  — — o n^  while  the 

e 

e”^  dependence  in  essentially  independent  of  the  value  of  n^.  Finally 


it  was  shown  that  resistivity  can  be  enhanced  due  to  the  density 

2 

fluctuation  with  ( lO”^  at  kX„  . 1 . A more  extensive 


treatment  of  this  subject  will  appear  elsewhere. 


6 


i . \ 


Acknowledgment 

We  would  like  to  express  our  thanks  to  Dr.  T.  Coffey 
for  pointing  out  the  data  and  the  possibility  of  inter- 
preting them  with  collisionless  processes.  Work  partially 
supported  by  NASA  Grant  and  by  ONR. 


References 


Amoldy,  R.  L,  and  Lawrence  W.  Choy,  Auroral  electrons  of  energy  less 
than  1 keV  observed  at  rocket  altitudes,  J.  Geophys.  Res.  JS, 

2187,  1973. 

Banks,  P.  H.,  C.  R.  Chappell,  and  A.  F.  Nagy,  A new  model  for  the 
interaction  of  auroral  electrons  with  the  atmosphere:  spectral 
degradation,  backscatter,  optical  omission,  and  Ionization, 

J.G.R.  1^59,  197^. 

Book,  D.  L.,  J.  P.  Boris,  and  K.  Haln,  Flux-Corrected  Transport  II: 

Generalizations  of  the  Method,  J.  of  Coop.  Phys.  3,  248,  197^. 

Berger,  M.  J.,  S.  M.  Seltzer  and  K.  Maeda,  Some  new  results  on  electron 
transport  In  the  atmosphere,  J.  Atmos.  Terr.  Phys.  591,  1974. 

DuBols,  D.  F.,  and  H.  V.  Goldman,  Parametrically  Excited  Plasma 
Fluctuations,  Phys.  Rev,  164,  207,  I967. 

Feldman,  P.  D.,  J.  P.  Doering,  and  J.  H.  Moore,  Rocket  measurements 
of  the  secondary  electron  spectrum  In  an  aurora,  J.  Geophys.  28, 
489,  1972. 

Feldman,  P.  D.,  and  J.  P.  Doering,  Auroral  electrons  and  the  optical 
emissions  of  nitrogen,  J.  Geophys.  Res.,  2d08,  1973* 

Matthews,  D.  L.,  M.  Pomgratz  and  K.  Papadopoulos,  Non-Linear 

production  of  suprathermal  tails  In  auroral  electrons,  J.  Geophys. 
Res.  123,  1976. 

Orena,  J.  H.,  J.  P.  Boris,  and  1.  Haber,  Optimization  techniques  for 
particle  codes.  Proceedings  Fourth  Conference  on  Numerical 
Simulation  of  Plasmas,  edited  by  J.  P.  Boris  and  R.  A.  Shanny 


11 


(U.S.  Government  Printing  Office,  Washington,  D.  C.,  1973)>  P-  526. 

Palmadesso,  P.  J.,  T.  P.  Coffey,  I.  Haber  and  K.  Papadopoulos,  Para- 
metric stabilization  of  relativistic  and  non-re lativistic  beam 
plasma  instabilities.  Bull.  Am.  Phys.  Soc.  IZJJ , 1975* 

Papadopoulos,  K.,  Nonlinear  stabilization  of  beam  plasma  interactions 
by  parametric  effects,  Phya.  Fluids,  1769>  1975* 

Papadopoulos,  K.  and  T.  Coffey,  Nonthermal  features  of  the  auroral 
plasma  due  to  precipitating  electrons,  J.  Geophys.  Res.  67k 
I97U. 

Papadopoulos,  K.,  and  T.  Coffey,  Anomalous  resistivity  in  the  auroral 
plasma,  J.  Geophys.  Res.  22^  J-558>  197^b. 

Papadopoulos,  K.,  M.  Goldstein,  and  R.  Smith,  Stabilization  of  electron 
streams  in  type  3 solar  bursts,  Astrophys.  J.,  190.  175»  197^* 

Rowland,  H.  L.,  and  K.  Papadopoulos,  Simulations  of  non-linear ly 
stabilized  beam  plasma  interactions,  to  be  published  (1977). 

Rowland,  H.  L.,  P.  J.  Palmadesso  and  K.  Papadopoulos,  One  dimensional 
d.c.  resistivity  due  to  strong  low  frequency  turbulence,  to  be 
published  (1977). 

Reasoner,  D L.,  and  C.  R.  Chappel,  Twin  pay  load  observations  of 
incident  and  backscattered  auroral  electrons,  J.  Geophys.  Res. 

18,  2176,  1975. 

Rees,  M.  H.,  and  K.  Maeda,  Auroral  electron  spectra,  J.  Geophys.  Res., 

28,  8391,  1973. 

Spiger,  R.  J.,  H.  R.  Anderson,  J.  S.  Murphree,  and  R.  F.  Loewenstein, 
Modulation  of  auroral  electrons  fluxes  in  the  frequency  range 
30  kHz  - 10  fflhz,  to  be  published,  1976. 


12 


Walt.  M.,  W.  M.  MacDonald,  and  W.  E.  Francis,  Penetration  of  auroral 
electrons  Into  the  atmosphere,  In  Physics  of  the  Magnetosphere, 
edited  by  R.  L.  Carovillano,  pp.  55^-555>  D*  Reldel,  Dordrecht, 
Netherlands,  1969* 


3 


ELECl 

< 

2 

f ; 

' I05 

L-L.l  J.U 

1 J 1.  1 J 1 LI 

I 2 3 4 6 8 10  20  30  40  60  80100 

ENERGY  (eV) 

Fig.  1 — Secondary  electron  spectrum  observed  by 
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